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The electronic structure of two V-based ladder compounds, the quarter-filled NaV^Os in the 
symmetric phase and the iso-structural half-filled CaV^Os, is investigated by ab initio calculations. 
Based on the bandstructure we determine the dielectric tensor e(u>) of these systems in a wide 
energy range. The frequencies and eigenvectors of the fully symmetric A g phonon modes and the 
corresponding electron-phonon and spin-phonon coupling parameters are also calculated from first- 
principles. We determine the Raman scattering intensities of the A 9 phonon modes as a function of 
polarization and frequency of the exciting light. All results, i.e. shape and magnitude of the dielec- 
tric function, phonon frequencies and Raman intensities show very good agreement with available 
experimental data. 

PACS numbers: 71.15.Mb, 63.20.Kr, 78.30.-j, 71.27.+a 



I. INTRODUCTION 

The vanadium-based ladder compounds NaV^Os and 
CaV^Os are interesting examples of systems where 
charge, spin, and lattice degrees of freedom are coupled 
to each other. Like materials with magnetically active Cu 
ions forming ladder-like structures they show unusual 
physical properties due to a strong interaction of all de- 
grees of freedom. The main building block of their unit 
cells is a ladder formed by V-0 rungs and V-0 legs, as it 
is shown in Fig. ^ Both materials crystallize in the or- 
thorhombic space group Pmmn (D^fJ with two formula 
units per unit cell. 




FIG. 1: The crystal structure of NaV 2 5 and CaV 2 5 . 
(a) the view along the ladders, (b) the view from the top. The 
dashed lines correspond to relatively weak bonds between the 
ladders. X indicates either Na or Ca, respectively. 



Electrons in these systems can move within a rung, be- 
tween rungs within a ladder, and between different lad- 
ders. The upper occupied electron states are formed by 
d xy electrons of V with a slight admixture of oxygen 2p 
orbitals. In NaV^Os, one d xy electron of V is shared by 
two sites within a rung which leads to a quarter-filled 
character— and makes this compound a dielectric, since 
the hopping of an electron between two rungs, which 
would produce a doubly occupied and an empty rung, re- 
quires more energy than can be provided by the hopping 
matrix element along the ladders. At a critical tempera- 
ture T c close to 35 K, NaV^Os undergoes a transition to 
a spin-gapped state, as first reported by Isobe and Ueda 
in Ref. 4j based on the analysis of magnetic susceptibility 
measurements. Evidence for such a transition was also 
obtained in Raman scattering experiments. 5 This behav- 
ior, in some aspects analogous to the spin-Peierls transi- 
tion observed in the inorganic chain compound CuGe03 
at Tsp — 13.5 K, 6 is accompanied by a disproportion 
of the V charges (4.5 ± S), a relatively large lattice dis- 
tortion, and the formation of an ordered charge pattern. 
At T C T c the ions are displaced from their equilibrium 
positions in the high-temperature phase by distances of 
the order of 0.05 A. These displacements which give ev- 
idence for a strong electron-lattice coupling have been 
observed in X-ray diffractionLSiS and can be estimated 
from infrared^ and Raman scattering experiments. 11 
They stabilize the zig-zag ordered phaseji2iiii being prob- 
ably the crucial element that determines the physics of 
the charge ordering in this compound. At the same 
time, the coupling to dynamical phonons induces strong 
charge fluctuations near T c which alter the spin-spin ex- 
change Jm^ Also the ordering is not completely static, 
even at low temperatures, since it is influenced by lat- 
tice vibrations, as was found in electron spin-resonance 
experiments. 15 At the same time, one could expect that 
charge ordering accompanied by a lattice distortion will 
show up, for example, in the dielectric function with de- 
creasing temperature. However, a comparison of optical 
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spectra taken at the low- and high-temperature phase, 
respectively, exhibits only small differences below T c ji& 
while Raman spectra reveal large changes in the elec- 
tronic background and show many new peaks which could 
have either magnetic or phononic origin when going be- 
low T c . n New peaks are also clearly seen in the low- 
temperature infrared spectral A full understanding of 
these new modes is still lacking. 

In CaV2C>5, each V ion has spin ^ and it can therefore 
be described by a generalized Heisenberg model of the 
spin-spin interaction. In this case, phonons directly in- 
fluence the exchange constants, and, therefore, lead to a 
modification of the magnetic properties. The spin gap in 
this compound, arising due to different exchange along 
the legs and along the rungs, is of the order of 500 K 
(0.05 eV)ii£ We mention that the Heisenberg spin-spin 
exchange parameters and, in turn, the strength of spin- 
phonon coupling depend on the electron on-site energies 
and the hopping matrix elements which form the band- 
structure and influence the dielectric function. For this 
reason, the experimental and theoretical investigation of 
the optical conductivity, the lattice dynamics, and the 
electron-phonon and spin-phonon coupling can provide a 
clue to the properties of NaV20s and CaA^Os and shed 
light on the origin of the phase transition in the former. 

NaV2C>5 and CaV205 have been the subject of in- 
tensive theoretical investigations, both by first principles 
based on density functional theory (DFT) and by model 
calculations ii&i2i22*2ii22i22i In first-principles calculations, 
the tight-binding linear muffin-tin orbitals (LMTO) 
method2i2^ in the atomic sphere approximation (ASA) 
and the linear combination of atomic orbitals^ (LCAO) 
approach, which directly included the Hubbard repulsion 
on the V sites, were applied to investigate the bandstruc- 
ture of NaV2C>5. More recently, Mazurenko et alSL com- 
bined DFT calculations with dynamical mean field theory 
aiming at understanding the insulating behavior of this 
compound in the high-temperature phase. The tight- 
binding LMTO method was also applied to CaV20s to 
obtain the spin exchange and effective four-band tight- 
binding model parameters for this compound^ In addi- 
tion, approaches based on quantum chemist Y y££&£±£L 
have been applied to NaA^Os. An interesting feature 
of the approaches used in Refs. OH] and 112 is that being 
based on the strong coupling of V and 01-p a orbitals, 
they attribute the phase transition to the ordering of the 
spin rather than the charge subsystem. 

While the bandstructures of NaV2 0s and CaV2 05 
are rather well understood on the first-principles basis, 
the analysis of their optical properties, lattice dynamics, 
electron-phonon coupling and Raman scattering still rely 
on various model assumptions. Having been well investi- 
gated experimentally, these properties require a detailed 
theoretical treatment which does not depend on such as- 
sumptions. Moreover, the outcome of DFT calculations 
is further used as realistic input parameters for model 
calculations such as Quantum Monte Carlo or Exact Di- 
agonalization techniquesii&M 



As far as NaX^Os is concerned, in this paper we con- 
centrate on the high-temperature phase where V ions 
in the rungs are equivalent. This enables us to under- 
stand its main properties and provides a starting point 
for the investigation of the low-temperature phase. Since 
in CaV20s no structural phase transition is observed, 
our treatment there holds at any temperature. The pa- 
per is organized as follows: In Section II, we describe the 
method of calculation and present the bandstructure for 
NaV2C>5 and CaV205 and related results like the den- 
sity of states and the charge density. We provide the cal- 
culated dielectric tensor components Su(uj) (i ~ x,y,z) 
for these two compounds in Section III. Section IV in- 
cludes theoretical phonon frequencies and eigenvectors, 
as well as electron-phonon and spin-phonon coupling pa- 
rameters. The phonon-induced changes in the dielectric 
function and the corresponding phonon Raman spectra 
will be presented in Section V. Finally, a summary of 
the results and suggestions for further investigations are 
given in the Conclusions. 



II. CALCULATIONS OF THE ELECTRONIC 
STRUCTURE 

A. Computational Methods 

All bandstructure calculations are performed within 
density functional theory (DFT) using the full-potential 
augmented planewaves + local orbitals (FP-APW+lo) 35 
formalism implemented in the WIEN2k codei2& Exchange 
and correlation terms are described within the gen- 
eralized gradient approximation (GGA)i21 The atomic 
sphere radii are chosen as 1.6 a.u. for V, 1.4 a.u. for the 
O atoms and Na, and 1.5 a.u. for Ca. In both compounds, 
all atomic positions have been relaxed starting from the 
experimentally measured ones as given in Ref. |3^ for 
CaV 2 5 and Ref. i for NaV 2 5 . In NaV 2 5 , for ex- 
ample, the shifts of the ions due to the relaxation of the 
structure are up to approximately 0.015 A (for 01 and 

02 in z direction) with an energy gain of roughly 30 meV 
per unit cell. 

Comparing the two materials, they have slightly dif- 
ferent lattice constants and ion coordinates. For exam- 
ple, since V is less positively charged in CaV205, the 
lengths of the V-0 bonds are slightly larger than those 
in NaV2 05. Specifically, the optimized values for the V- 

03 and V-Ol bond lengths are 1.62 (1.67) A and 1.82 
(1.85) A in NaV 2 5 (CaV 2 5 ). At the same time, the 
spacing between Ca and the 01 plane (2.39 A) is smaller 
than the Na-Ol plane distance in NaV205 (2.44 A), since 
the bigger Ca ion exhibits a stronger Coulomb interaction 
with oxygen compared to Na. The geometry relaxation 
allows to make the calculations not directly relying on 
the experimentally measured structural data and thereby 
leads to small quantitative differences compared to the 
NaV2 05 calculations performed by Smolinski et al& At 
the same time, this provides the energy scale related to 
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unit cell distortions as it is realized for example in the 
low-temperature phase. 

We do not include correlational effects by using an 
LDA+f/ approach in our calculations, but we have es- 
timated the Hubbard U for NaV20s and CaV2C>5 from 
our data by the following procedure. Similarly to what 
is described in Ref. |3J, we have added a small amount 
of electronic charge to the system and estimated U from 
the resulting change of the V bands. Charge neutral- 
ity was accounted for by two different procedures: When 
putting the positive charge on the Na sites, the resulting 
U was estimated to be 2.8 eV for NaV205 (averaged over 
the Brillouin zone (BZ)) in good agreement with Smolin- 
ski et al. 3 We preferred, however, to provide the positive 
charge in terms of a uniform background, which leads to 
a U value of 2.45 eV for both compounds, demonstrating 
that U only weakly depends on the ion's surrounding. 
This procedure has the advantage, that the band energy 
shifts are much more uniform with respect to different k 
points of the BZ than when the additional positive charge 
is located at the Na sites (differences of hundredths of an 
eV in the former case compared to tenths of an eV in the 
latter case). 



B. Bandstructure and density of states 

The bandstructures and densities of states of NaV2C>5 
and CaV2C>5 are shown in Figs. and |3J respectively. 
Focusing on NaV20s first, the bands 3-4 eV below the 
Fermi level are due to 0-2p states. Their smallest dis- 
tance to the valence band minimum (at the T point) will 
be denoted as E g later in the text. The unoccupied states 
exhibit mainly V character with a small admixture of 
oxygen. The bands intersecting the Fermi level in k y 
direction are formed by the bonding combination of V 
orbitals. Their dispersion is due to hopping along the 
ladders while the splitting of this pair of bands at the 
r point is due to inter-ladder hopping. The two bands 
just above the Fermi level originate from the antibonding 
combination of V-d xy states. The situation for CaV205 
is similar, where some quantitative differences will be dis- 
cussed below and in the context of its optical properties. 

The bands can be mapped onto a tight-binding model 
with the one-ladder parameters t± of the in-rung hop- 
ping, and tii representing the hopping along the lad- 
der. The theoretically determined values for NaV2C>5 
are t± — 0.387 eV and tn = 0.175 eV, which is close to 
the data of Ref.0 For CaV 2 5 , we obtained in the same 
way t± = 0.321 eV and tn — 0.143 eV, in agreement with 
the results of Korotin et alM^ who applied the LDA+C7 
technique in their calculations. Compared to their re- 
sults, our hopping matrix elements are slightly increased, 
which is due to the changes in the interatomic distances 
as a result of the structural relaxation. We note that both 
t± and tn, are smaller for CaV205 than for NaV205. At 
the same time, the splitting of the bonding bands in the 
r point arising from the inter-ladder hopping is much 
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FIG. 2: Bandstructure and density of states (in states per 
unit cell and eV) of NaV 2 5 . 



larger in CaV 2 5 (0.49 eV) than in NaV 2 5 (0.23 eV) 
since the distance between the ladders is smaller in the 
former. The corresponding inter-ladder hopping matrix 
elements between the closest V atoms of neighboring lad- 
ders, ti, are 0.13 and 0.24 eV for NaV 2 5 and CaV 2 5 , 
respectively. When the lattice is deformed by a displace- 
ment of ions corresponding to a phonon mode, the tight- 
binding parameters as well as the on-site energies change. 
This kind of electron-phonon coupling will be discussed 
below. 

The influence of the Hubbard term U on the properties 
of V-based ladder compounds is widely discussed in the 
literature. It is important to mention that the enhanced 
electron correlation when accouted for by the Hubbard 
parameter reproduces the semiconducting behavior with 
the charge gap close to 2|tj_| , 18 Correspondingly, the 
dispersion along the y-axis of the band derived from the 
bonding combination of V-d xy orbitals will be n rather 
than 2-7T periodici^MS At the same time, we shall see be- 
low that the physical properties determined by the elec- 
tron density, are not strongly influenced by the Hubbard 
repulsion and can be described reliably within DFT. 

To illustrate the charge density distribution within the 
unit cell, we exemplarily present the electron charge den- 
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FIG. 3: Bandstructure and density of states (in states per 
unit cell and eV) of CaV2 0s. 

sity for NaA^Os in Fig. 0] in two perpendicular planes. 
The upper panel clearly shows the role of the unit cell 
asymmetry on the charge density arising due to the pres- 
ence of the apical oxygen 03. This asymmetry, on the one 
hand, leads to a strong Holstein-like electron-phonon cou- 
pling, and, on the other hand, diminishes the overlap of 
the V orbitals with 01 and 02 states, thus decreasing the 
hopping matrix elements and correspondingly the com- 
ponents of the dielectric tensor. In addition, the lower 
panel shows the preferred orientation of the in-ladder 
oxygen states and a relatively small overlap of the or- 
bitals from different ladders. 



III. DIELECTRIC FUNCTION 

Figures |S] and [S] present the real and imaginary parts of 
the diagonal dielectric tensor components, Re £u(uj) and 
Im Ejj(a;), in the experimentally measured range, where 
the Cartesian index i corresponds to the light polariza- 
tion directions. Im Su(lo) was calculated within the Ran- 
dom Phase Approximation (RPA), based on the Kohn- 
Sham orbitals including a life-time broadening of the op- 
tical inter-band transitions of 0.1 eV, while Re sa(u)) is 
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FIG. 4: Valence electron density of NaV20s in the (a;, z) plane 
(upper panel) and the (x,y) plane (lower panel), which both 
contain the V positions. The logarithmic contour lines range 
from 0.07 e/A 3 to 2.3 e/A 3 . 



obtained by Kramers-Kronig transformation. 

Let us first discuss the in-plane response of NaV205. 
The most interesting feature of the xx component (light 
polarized along the rungs) is a strong peak at 1.03 eV 
in agreement with experiment ^S^LiS An analysis of the 
interband momentum matrix elements at different elec- 
tron wavevectors k v shows that the first peak in the xx 
response arises due to transitions between the bonding 
and antibonding band states within one V rung. The en- 
ergy of the peak is larger than 2t± because of the band 
dispersion along the y axis and can be estimated to be 
2(i_i_+ii|), which is 1.12 eV. Since the transitions between 
d xy orbitals have very small matrix elements due to the 
large V-V distance, the admixture of 0\-p y states in the 
antibonding state is responsible for a sizable intrarung 
transition matrix element. We mention here that this 
admixture provides some support to the arguments of 
Refs. l3ll andl32l It rapidly decreases, however, with the 
increase of k y due to the corresponding decrease of the 
oxygen contribution, which was first noticed in Ref. U 
For this reason the peak mostly originates from transi- 
tions in the vicinity of the T point. 

The yy component (light polarized along the legs) is 
dominated by a double-peak structure at 1.27 eV and 
1.78 eV, respectively. It is considerably weaker than the 
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FIG. 5: Diagonal components of the dielectric function eu(u>) 
ofNaV 2 5 . 
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FIG. 6: Diagonal components of the dielectric function eu(uo) 
ofCaV 2 5 . 



xx response. The analysis of the bandstructure reveals 
that the shoulder at 1.3 eV comes from in-rung transi- 
tions. These can contribute to s yy (u) since at finite k x 
values the in-rung states are neither odd nor even with 
respect to the x — > —x transformation, and, therefore, 
can couple to light polarized along the y-axis. The broad 
maximum at u) = 1.8 eV is related to transitions from 
the bonding V-d xy states to 02-p x states admixed to V- 
d^-derived orbitals at approximately 1.7 eV above the 
Fermi level. 

In both polarizations, the peaks occurring at higher 
energies {w > 3 eV) originate from transitions between 
0-2p and V orbitals. For example, the broad feature 
around 3 eV in the xx spectra arises from transitions be- 
tween 01-2p z states at -3.3 eV and bonding V-d xy states 
around Ep. Comparing to experiment, we want to point 
out that all theoretically obtained features reproduce the 
corresponding experimental findings very well^& 

In CaV2 05, the character of the transitions is, in gen- 
eral, the same as in NaV2 0s. The fact that Ca provides 
one more valence electron to the system compared to 
Na has two obvious effects: First, the optical response is 
stronger compared to NaV2 0s. Second, the peak present 
in NaA^Os around 3 eV is missing, since the bonding V- 
d xy states are occupied and thus don't provide final states 
for the transitions. At the same time, since the in-ladder 
hopping matrix elements are smaller in CaV2 05 than 
in NaV205, the spectrum of unoccupied band states in 



CaV205 is denser, as can be seen in Figs. [21 and [31 For 
this reason, the low-energy part of the dielectric func- 
tion (u < 1 eV) of CaV2C-5 shows a more complicated lj 
dependence than that for NaV20s since more interband 
transitions are allowed at this spectral range. A compar- 
ison with experiment for this compound which can also 
provide an experimental test of the applicability of the 
DFT for the description of optical properties of half-filled 
ladder compounds is not possible at present since to the 
best of our knowledge no measured data of the dielectric 
function of CaV20s single crystals are available. 



IV. LATTICE DYNAMICS 

A. Phonon modes 

For the calculation of the A g phonon modes we ap- 
plied the frozen-phonon approximation. To this extent, 
the atomic positions have been moved from their equilib- 
rium. Four displacements (two in positive and negative 
direction, respectively) for each degree of freedom have 
been taken into account. The resulting forces were used 
to obtain the energy hyper-surface according to the pro- 
cedure described in Ref. ^3 and to set up the dynamical 
matrix. Since the harmonic fully-symmetric ion vibra- 
tions do not change the occupancy of the V sites from 
single (of fractional) to double, the electron correlations 
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V z 


Na z 


Olz 


02 x 


02 z 


03x 


03 z 


Refs. 44.46 


This work 


970 


969 


996 


0.04 


0.25 


-0.01 


0.01 


0.01 


-0.00 


-0.05 


-0.43 


V-03 stretching 


V-03 stretching 


con 
ooU 


0o4 


512 


0.16 


-0.07 


0.01 


0.16 


0.45 


0.06 


-0.05 


-0.01 


V-02 stretching 


V-02 stretching 


450 


448 


467 


0.42 


-0.06 


0.01 


0.15 


-0.21 


-0.05 


-0.11 


0.01 


V-Ol-V bending 


V-Ol-V bending 


422 


423 


414 


-0.19 


-0.17 


0.00 


0.44 


-0.08 


-0.18 


0.18 


-0.14 


03-V-02 bending 


01 z +03-V-02 bend. 


304 


304 


308 


0.02 


0.19 


-0.04 


0.18 


-0.07 


0.39 


0.19 


0.09 


(J3-V-U2 bending 


U3-V-U2 bending 


230 


233 


232 


0.17 


0.06 


0.03 


-0.13 


0.06 


-0.16 


0.42 


-0.00 


03-V-02 bending 


03- V-02 bending 


178 


179 


176 


-0.02 


0.28 


-0.40 


0.10 


0.04 


-0.22 


-0.06 


0.18 


Na || c 


Na || c 


90 


90 


111 


-0.04 


0.30 


0.42 


0.13 


0.03 


-0.16 


-0.08 


0.18 


chain rot. 


chain rot. 



TABLE I: Calculated frequencies and eigenvectors of the A 9 phonon modes of NaV20s compared to experiment. 
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0.06 


0.24 


0.01 


0.02 


0.02 
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-0.42 


V-03 stretching 


V-03 stretching 
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539 
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-0.01 
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-0.01 


0.11 


0.49 


0.07 
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0.00 


V-02 stretching 
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0.39 


-0.09 


-0.05 


-0.08 


-0.00 
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421 
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-0.34 


-0.12 


0.02 


0.27 


-0.05 


-0.17 


0.18 


-0.13 


03-V-02 bending 01 z +03-V-02 bending 


337 


? 
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0.20 


0.07 


0.10 
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0.20 
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0.05 
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265 


0.19 
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0.09 


-0.31 


0.25 
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03- V-02 bending 


03-V-02 bending 
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138.6 


201 


0.02 


0.21 


0.39 


0.03 


0.06 


-0.31 


-0.05 


0.17 


Ca || c 


chain rot. 


91 


90 


106 


-0.05 


0.32 


-0.39 


0.15 


0.01 


-0.15 


-0.06 


0.19 


chain rot. 


Ca | c 



TABLE II: Calculated frequencies and eigenvectors of the A 9 phonon modes of CaV20s compared to experiment. 



do not significantly influence the elements of the dynam- 
ical matrix. At the same time, the correlation effects 
could be more important for the anharmonic terms re- 
lating the lattice forces and ion displacements. 

The frequencies of the A g modes for NaV2 0s and 
CaV2C>5 are presented in Tables [I] and |nl respectively, 
and compared to experimental data. The corresponding 
eigenvectors presented in these Tables, e a ^, are related 
to the real displacements u" by 

e Q c = const, x u"\/M Q , (1) 
and are normalized as 

N 

E e «C = l> (2) 

Q = l 

where a enumerates the ions with mass M a , N = 16 
is the number of ions per unit cell, and £ indicates the 
phonon mode. In order to visualize the lattice distortions 
according to the phonon eigenvectors, the corresponding 
atomic displacements of NaV2C>5 are shown in Fig. 
where the eigenvector components of Table [|] refer to the 
equivalent positions labeled in the figure. The phonon 
eigenvectors of CaV2C>5 are very similar, up to some dif- 
ferences discussed later. 



For the eigenfrequencies of the NaV2 0s A g modes, 
good agreement with experiments^*^^ is found, with 
deviations smaller than 5%. Only for the lowest-energy 
mode the difference is larger. In full agreement with 
the results of Refs. HIE and the eigenvector of the 
996 cm" 1 mode represents a stretching between V and 
the apical oxygen. Also the V-02 stretching of the 
512 cm -1 mode and the V-Ol-V bending of the 467 cm -1 
mode as suggested in Ref. ^| are verified, where we, how- 
ever, find admixtures of an 01 z motion for the former 
and an 02 x movement for the latter, respectively. In the 
eigenvector of the 414 cm- 1 mode, the z-displacement of 
the in- rung oxygen 01 is dominating, while this mode is 
described as pure 03- V-02 bending in Ref. 0. At the 
same time, our theoretical frequencies are much closer to 
experiment (2% deviation) than the calculated frequen- 
cies in Ref. S (9% difference). For the 308 cm" 1 , the 
232 cm -1 , the 176 cm.- 1 and the 111 cm -1 mode, the 
agreement of our results with the assignment of Popovic 
et a/4& is good. However, in most of the modes we find 
a more pronounced involvement of 01 compared to the 
interpretation of experimental results. 

For the eigenfrequencies of the CaV205 A g modes, the 
agreement between theory and experiment^ is very good 
for the modes above 400 cm -1 . For these vibrations, the 
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996 cm 1 : V-03 stretching 
x 




467 cm'' : V-Ol-V bending 



414 cm" 1 : 01 z + 03-V-02 bending * 03 



996 cm 1 : V-03 stretching 
x 




512 cm" 1 : V-02 stretching 



467 cni - ' : V-Ol-V bending 



414 cm" 1 : 01 z + 03-V-02 bending ^<j3 ^ 




FIG. 7: Eigenvectors of NaV 2 5 . 

experimental assignment with respect to their symmetry 
is unambiguous. The lower frequency of the apical oxy- 
gen vibration in CaV2C>5 compared to NaV20s is due 
to larger interionic distances, and hence, smaller force 
constants. A change of 5% can be estimated within the 
Coulomb picture from the different Vanadium charges in 
NaV2C>5 and CaA^Os (i.e. 4.5 and 4, respectively), which 
is in qualitative agreement with experiment. 



The two lowest-frequency modes are swapped when 
Na is replaced by Ca, i.e. the in-phase motion of Ca 
with the other atoms of the ladder (chain rotation) has 
higher energy than the out-of-phase vibration, where Ca 
vibrates || c in opposite direction to the ladder. The fre- 
quency of the chain rotation mode (201 cm -1 ) is roughly 
twice as high as in NaA^Os, representing a difficulty in 
the interpretation of the measured modes of CaV205 
which in Refs. and was done in comparison with 
the phonons of NaV20s. In Ref. the 138 cm -1 mode 
of the unpolarized spectrum was interpreted as A g vi- 
bration since its frequency compared to that of the (Na 
|| c) mode scales as the inverse square root of the cor- 
responding masses. From our analysis we conclude, that 
the measured 138 cm -1 mode is not an A g vibration. The 
physical origin of the frequency shift of the chain rotation 
mode is related to the stronger inter-ladder interaction 
in CaV2C>5 which is also reflected in the enhanced tight 
binding parameter t, (see Section IV B). On the other 
hand, the higher Ca mass is mainly responsible for the 
decrease of the "Ca || c" mode frequency. It is still un- 
clear, however, whether the calculated 201 cm -1 mode 
corresponds to the experimentally observed 213 cm -1 
vibrationaa or, more probable, to a mode at 235 cm -1 
(Ref. l4r?) (or 238 cm -1 according to Ref. i45h which has 
been assigned as an A g mode in both papers. In the 
latter case our calculated frequency of 265 cm -1 can be 
related to the measured 282 cm -1 vibration^ A similar 
problem concerns the interpretation of one more A g mode 
somewhat above 300 cm -1 , where no clear experimental 
assignment is available in Ref. In this context, the 
interpretation of the Raman scattering intensities could 
be helpful for an unambiguous assignment as will be dis- 
cussed in the next section. 



B. Electron-phonon and spin-phonon coupling 

When ions are shifted from their equilibrium posi- 
tions, the changes of the bandstructure are a measure 
for the electron-phonon interaction. Two types of cou- 
pling can be considered: The first one is the Holstein 
coupling, where the site energies change with the ion dis- 
placements, while the other one is due to changes of the 
hopping parameters. At the same time, the exchange 
path, which is formed by transitions between different 
sites, is also influenced by the phonons. This effect re- 
sults in spin-phonon coupling. To investigate the type 
of the electron-phonon coupling and its strength for each 
mode, in Tables IIIII and IIVI we display the changes of 
various model parameters with the corresponding ion dis- 
placements u" (up to ~ 0.05 A). These are the hopping 
parameters tu, t±, and ti, the charge transfer gap E gi 
and the exchange parameters J\\ (along the ladders) and 
Ji (within the rungs). Jy and Jj_ can be estimated as 
~ i| /Eg and ~ t\/E g , respectively^ In this context we 
introduce the dimensionless phonon coordinate Q by the 
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relation 



M a U>£ 



(3) 



It turns out that only the higher-frequency phonons 
considerably modulate the one-ladder parameters t± and 
t||. The changes of the hopping matrix elements corre- 
sponding to Q = 1 are below 0.025 eV, while they can be 
much larger in the energy shifts 5E g . Therefore the main 
mechanism of electron-phonon coupling can be assigned 
to Holstein-like interaction. Our results also allow to esti- 
mate the strength of spin-phonon coupling arising due to 
the phonon-induced modulation of the exchange parame- 
ters. The corresponding relative changes are summarized 
in Tables HTT1 and llVl 

The changes in the matrix elements can be under- 
stood as a result of the altered interionic distances and 
the electronic on-site energies. For example, the large 
phonon-induced decrease of t\\ in the 512 cm" 1 vibration 
of NaV2 05 is a consequence of the larger V-0 distance 
in the same leg (see Fig. 01 . At the same time, the V-0 
distance along the x axis changes from 0.274 A to 0.309 A 
at Q = 1. The decrease of t± for the 308 cm -1 mode is 
due to an enhanced z-axis distance between V and the 
in-rung oxygen 01 by 0.02 A going from the relaxed to 
the distorted structure with Q = 1. As a consequence, 
also the energy difference of the V and 01 orbitals is 
increased, diminishing the hopping parameter t±. The 
very strong modulation of the inter-ladder hopping t i by 
the 512 cm -1 mode is due to a zig-zag like deformation of 
the legs in the (x, y) plane, i.e. a vibration of neighbor- 
ing vanadium and oxygen atoms of one leg in opposite 
direction, and hence an increase of this V-O2 distance. 

As can be seen in Tables Hill and llVl the biggest change 
for both in and t±, in NaV205 is caused by the V-02 
stretching mode, while in CaV20s this mode leads to 
the most pronounced change only in tii , but the largest 
modulation of t± is due to the V-Ol-V mode (446 cm" 1 ). 
The reason for this can be found in the displacement of 
the in- rung oxygen 01 along the z axis, which is much 
larger compared to NaV2 05. For both compounds, the 
biggest effect on E g is observed for the V-03 stretching 
mode. 

We emphasize here that both, electron-phonon cou- 
pling (leading to a modulation of E g and the hopping 
matrix elements) and spin-phonon coupling (leading to a 
modulation of J||) in NaV205 and CaV20s, are consid- 
erably strong^. For this reason the lattice distortion in 
the low-temperature phase of NaV205 can be related to 
total-energy changes originating from charge as well as 



,31.32 



degrees of freedom. 



V. RAMAN SCATTERING 

With the knowledge of the phonon modes and the di- 
electric functions we can calculate the phonon Raman 



spectra of NaV205 and CaV205. For this purpose we 
use the approach developed in Ref. l43l where at a given 
exciting light frequency luj the total Raman intensity In 
at temperature T in arbitrary units is: 

= V(n s ( Wc ) + 1) (1|_JLQ|0> L(^,^ c ,r) 



C 



dQ 



(4) 

Here lur is the Raman shift, the Cartesian indices ii 
correspond to the polarizations of incident and scat- 
tered light, which are the same due to the orthorhombic 
symmetry of the crystal. |1) and |0) denote the one- 
phonon and the phonon-less states, respectively, and Q 
is the operator of the phonon coordinate, nsj^) = 
1 / (exp (huj( /T) — 1) is the phonon Bose distribution func- 
tion, and L(ujr,ijJq,T) is the Lorentzian shape of the 
phonon line with a broadening T, which was chosen to 
be 25 cm -1 for all modes. 

The total Raman intensity, i.e. the sum over all phonon 
contributions, is presented in Figs. |51 and 03 for Nav^Os 
and CaV205, respectively, for an incident light energy of 
2.41 eV (A = 514.5 nm), which is used in the Raman ex- 
periments available in the literature. In the xx polariza- 
tion seven out of eight modes are clearly visible, only the 
232 cm -1 vibration has negligible intensity, in excellent 
agreement with experimentsiAi*24 We note that also the 
relative peak heights are fully reproduced. This scatter- 
ing geometry exhibits the highest intensity for all modes 
except the highest one, which dominates the zz polarized 
spectra. The intensity of all other modes in this polariza- 
tion is two orders of magnitude smaller and hence hardly 
visible in the measured spectra. The only exception is 
the 110 cm" 1 mode, where theory cannot reproduce the 
experimentally observed sharp peak. One possible expla- 
nation could be provided by the extremely pronounced 
resonance behavior of most of the vibrations as will be 
discussed below. The yy polarization exhibits intensities 
in between the magnitudes of the xx and zz counterparts. 
In accordance with experiment the 467 cm -1 vibration 
is absent in this scattering geometry. 

For CaV205 the situation is similar. The highest en- 
ergy mode exhibits the highest intensity in zz polariza- 
tion. Between 400 and 600 cm -1 the xx intensities are 
dominating. Only below 400 cm -1 , the yy spectra are 
comparable in magnitude or even bigger. 

The Raman scattering intensity is governed by the 
change of the crystal polarizability with the nuclei vi- 
brating around their equilibrium positions. The depen- 
dence of the dielectric function on the phonon coordinate 
arises due to two main reasons, i.e. the Q dependence of 
the momentum matrix element and the Q-dependence of 
the interband transition energy^ The latter contribu- 
tion leads to a stronger resonance behavior of the Ra- 
man intensity than the former. To illustrate the influ- 
ence of the lattice vibrations on the dielectric function, 
in Fig.^Jwe present \deu(ui) /dQ\ 2 for the phonon modes 
of NaV205. In the energy range ui > 2.5 eV the deriva- 
tive of the xx component dominates over \de yy {ijj) / dQ\ 2 
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oJth (cm 


- 1 ) 


8t\\/5Q (eV) 


8t ± /8Q (eV) 


5U/8Q (eV) 


SE g /8Q (eV) 


(SJ/J^/SQ 


(8J/J)±/8Q 


996 




-0.0021 


-0.0099 


0.0056 


-0.0936 


0.0190 


0.0449 


512 




-0.0210 


0.0213 


-0.0439 


0.0905 


0.2952 


-0.0773 


467 




-0.0016 


-0.0115 


0.0109 


-0.0037 


0.0052 


0.0478 


414 




-0.0110 


-0.0069 


-0.0142 


-0.0010 


0.1284 


0.0294 


308 




0.0000 


-0.0167 


-0.0253 


0.0554 


0.0003 


0.0967 


232 




-0.0012 


-0.0012 


0.0257 


-0.0340 


0.0085 


0.0013 


176 




0.0019 


0.0061 


0.0193 


-0.0185 


-0.0130 


-0.0225 


111 




0.0008 


0.0028 


0.0132 


-0.0117 


-0.0055 


-0.0107 




undistorted: tii 


= 0.175 eV, t± = 0.387 eV, U = 0.117 eV, E a = 2.565 eV 



TABLE III: Parameters of electron-phonon and spin-phonon coupling for the A 9 eigenmodes of NaV20s. 



a>th (cm -1 ) 


8t\\/8Q (eV) 


St ± /8Q (eV) 


8U/8Q (eV) 


8EJ8Q (eV) 


(SJ/J^/SQ 


(SJ/J) ± /SQ 


900 


-0.0029 


-0.0089 


0.0078 


-0.1286 


0.0316 


0.0470 


516 


-0.0113 


0.0120 


-0.0253 


0.0503 


0.1863 


-0.0565 


446 


-0.0013 


-0.0232 


0.0139 


-0.0228 


0.0052 


0.1399 


412 


0.0056 


-0.0094 


0.0223 


-0.0239 


-0.0815 


0.0540 


307 


0.0016 


-0.0017 


-0.0035 


0.0054 


-0.0141 


0.0185 


265 


0.0071 


-0.0050 


-0.0337 


-0.0048 


-0.0895 


0.0285 


201 


-0.0019 


0.0035 


0.0278 


-0.0080 


0.0269 


-0.0217 


106 


0.0013 


-0.0055 


-0.0110 


-0.0009 


-0.0227 


0.0317 


undistorted: ty 


=0.143 eV, tj 


= 0.321 eV, U = 0.244 eV, E a = 2.882 eV 



TABLE IV: Parameters of electron-phonon and spin-phonon coupling for the A 9 eigenmodes of CaV^Os. 



for all modes. In the low energy region the latter shows 
a dramatic resonance behavior for the 232 cm -1 and the 
308 cm -1 vibration, and is comparable or even slightly 
larger in magnitude than \de xx (uj)/dQ\ 2 around 2 eV for 
the next three modes higher in energy. All modes but 
the highest one exhibit extremely strong resonances in 
the infrared in either of the two polarizations. For the 
highest mode, the zz polarization is dominating, while 
only below 2 eV the yy component becomes more pro- 
nounced. We thus predict much higher Raman scattering 
intensities for incident light frequencies both higher and 
lower than those used in experiments so far. We mention 
in this context, that Fischer et ali* observed strong rel- 
ative intensity changes in xx geometry when measuring 
the low temperature phase using different photon ener- 
gies. 



A similarly strong resonance behaviour is also expected 
for CaV205. While with the exception of the high- 
est mode the intensities of the zz component are small 
and the influence of the photon energy is negligible in 
NaV2C>5, this scattering geometry is slightly richer for 
CaA^Os, with generally slightly higher relative intensi- 
ties and an even more pronounced resonant behavior of 
the highest mode. 



VI. CONCLUSIONS 



In this paper, we have studied optical properties and 
lattice dynamics of NaV20s and CaV20s in the Pmmn 
phase from first principles. The calculations are based 
on the theoretically optimized crystal structures obtained 
within the generalized gradient approximation. Effective 
bandstructure parameters have been extracted by map- 
ping our results onto a tight-binding model. We have 
obtained the Hubbard repulsion U on the V sites being 
approximately 2.45 eV, both for NaV20s and CaV205. 
The dielectric functions have been determined within the 
random phase approximation and are in very good agree- 
ment with available experiments. Our results show that 
the 1 eV peak in the xx component arises due to tran- 
sitions between the bonding and the antibonding combi- 
nation of V d xy orbitals within one rung. By diagonaliz- 
ing the dynamical matrix we have obtained the phonon 
frequencies for the fully symmetric vibrations which are 
in very good agreement with measured data. With the 
knowledge of the phonon eigenvectors and the changes of 
the bandstructure caused by the phonon modes we have 
estimated the parameters of electron-phonon and spin- 
phonon coupling for both compounds. We find that the 
strongest contribution to the electron-phonon coupling 
comes from the phonon-modulation of the charge trans- 
fer gap Eg. At the same time, other effects arising from 
altered hopping matrix elements can be important. Fi- 
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FIG. 8: Raman intensity Ir for NaV^Os in three different 
geometries at u = 2.41 eV (A = 514.5 nm) and T = 300 K. 



nally, we have calculated the phonon Raman spectra of 
these compounds and analyzed the frequency-dependent 
dielectric function modulated by the ion displacements 
according to the lattice vibrations. On this basis we pre- 
dict a strong resonance behavior for both NaA^Os and 
CaV 2 5 . 



As an outlook for further investigations, this detailed 
analysis provides a basis for comparison with the low- 
temperature phase. Moreover, the Hubbard parameters 
extracted in this work can be used as input for further 
calculations by, e.g., Exact DiagonalizationS&SH or Quan- 
tum Monte Carlo simulation^ 
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